Two-channel model description of confinement-induced Feshbach molecules 
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Using a two-channel model, we investigate theoretically the binding energy of confinement-induced 
Feshbach molecules in two- and one-dimensional ultracold atomic systems, near a Feshbach reso- 
nance. We show that the two-channel prediction will evidently deviate from the simple single-channel 
theory as the width of Feshbach resonances decreases. For one-dimensional system, we perform a 
full two-channel calculation, with the inclusion of bare interatomic interactions in the open channel. 
Away from the resonance, we find a sizable correction to the binding energy, if we neglect incorrectly 
the bare interatomic interactions as in the previous work [Dickerscheid and Stoof, Phys. Rev. A 
72, 053625 (2005)]. We compare our theoretical results with existing experimental data and present 
predictions for narrow Feshbach resonances that could be tested in future experiments. 



I. INTRODUCTION 

The confinement-induced resonance (CIR) is one of the 
most intriguing phenomena in low-dimensional ultracold 
atomic systems, and has attracted a great deal of interest. 
It was first predicted by Olshanii using a single-channel 
model in 1998 when considering two-body collisions in 
one-dimensional (ID) harmonic waveguides [lj. Later, 
this study was extended to two-dimensions by Petrov and 
co-workers in 2000 To date, there are a number 

of experimental confirmations of the existence of CIR in 
both ID and 2D setups. 

To confirm experimentally CIR, it is convenient to use 
a Feshbach resonance to tune the interatomic interactions 
d, [f| and to measure the binding energy of the resulting 
confinement-induced Feshbach molecules (CIFMs). The 
existence of CIFMs in a ID Fermi gas of 40 K atoms was 
first detected by Moritz et. al. Using CIR in a 

ID Bose gas of 133 Cs atoms @|, Haller et. al. realized 
the crossover from the Tonks-Girardeau [H, HJ to super- 
Tonks- Girardeau regime [Icj . Interestingly, in such ID 
experiments [ll| , an anomalous splitting of the CIRs ap- 
peared when an anisotropy in the transverse confinement 
was introduced, which could not be explained by general- 
izing Olshanii's theory to anisotropic transverse confine- 
ment [12]. But soon, it was realized that the splitting 
of the CIRs was resulted from the coupling between the 
center-of-mass (COM) and relative motions of two incom- 
ing atoms due to the anharmonicity of the trap [HI, [l4j . 
Besides these ID experiments, CIR was also recently ob- 
served in a 2D Fermi gas of 40 K atoms [l5j. The in- 
teraction between atoms was tuned by using a magnetic 
Feshbach resonance at Bo = 224. 2G, and the binding en- 
ergy of the CIFMs was measured. Compared with the 
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prediction of the single-channel theory 0, however, 
the observed binding energy is larger than the theoretical 
prediction by 4 kHz in magnitude fl5] . This discrepancy 
may be understood from the picture of fermionic polaron, 
as suggested by Schmidt et al. [16|. Alternatively, the 
discrepancy may also be solved by introducing an energy- 
dependent interaction, as shown experimentally for 40 K 
atoms near another Feshbach resonance Bq — 202. 1G 

Here, we aim to study CIR and CIFMs by using a 
two-channel model, with the inclusion of the effect of 
bound molecule states in the closed channel. This is a 
physically more realistic description for Feshbach reso- 
nances and therefore should provide a better description 
for CIFMs which are measured experimentally. It is well- 
known that the two-channel model yields the same the- 
oretical predictions as the single-channel model in the 
limit of broad Feshbach resonances (l9l . [20| . The exper- 
imentally utilized resonances for 40 K and 133 Cs atoms 
are broad. As a result, our two-channel model descrip- 
tion may not give improved understanding of the existing 
CIR measurements. However, our two-channel results 
should provide a useful guide for future CIR experiments 
on relatively narrow Feshbach resonances. 

We note that, a two-channel calculation was previously 
carried out for ID CIFMs by Dickerscheid and Stoof in 
2005 [18]. However, in their study the bare interatomic 
interactions were neglected. It was shown then the two- 
channel model gives a better agreement with the ID ex- 
perimental data than the single-channel model |18| . 

In this work, we present a two-channel calculation for 
the binding energy of CIFMs in 2D systems under an 
axially harmonic confinement. Near the Feshbach res- 
onance, the channel coupling will dominate the contri- 
bution to the binding energy of the dressed molecules, 
which allows us to neglect the bare interatomic inter- 
actions in the open channel. At this resonance limit, we 
find that the two-channel result coincides with the single- 
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channel theory for a broad Feshbach resonance, but will 
evidently deviates from the single-channel theory as the 
resonance width decreases. We compare our theoreti- 
cal calculation with the recent 2D 40 K experiment [l7j . 
in which a broad Feshbach resonance at Bo = 202. 1G is 
used to control the interatomic interactions. We find that 
the two-channel theory is in agreement with the single- 
channel result, as expected, and both agree well with the 
experimental data at the resonance. In order to demon- 
strate the in-equivalence between the two-channel and 
single-channel theories, we predict the binding energy of 
2D CIFMs for two relatively narrow resonances, 23 Na at 
B a = 907G and 87 Rb at B a = f007.4G. 

We also perform a full two-channel calculation for ID 
CIFMs. We find that the bare interatomic interactions 
will give rise to a sizable correction to the binding energy 
away from Feshbach resonances. This bare interatomic 
interaction was neglected in the previous treatment [l8| . 
By restoring the bare interatomic interactions, the two- 
channel result coincides with the single-channel predic- 
tion, as we anticipate for a broad resonance. Thus, we 
realize that the previous better agreement claimed by 
Dickerscheid and Stoof is not convincing. 

The paper is arranged as follows. We first present the 
two-channel Hamiltonian in Sec. |TTJ Then the ansatz of 
two-body wavefunctions for 2D and ID systems is con- 
structed in Sec. Mil In Sec. IIV1 we solve the two-body 
problems and calculate the binding energy of CIFMs. In 
Sec. El we report the binding energy as a function of 
the width of Feshbach resonances and show how the two- 
channel result deviates from the single-channel predic- 
tion as the resonance width decreases. We also compare 
our theoretical results with the recent experiments and 
discuss in detail the correction of the bare interatomic 
interactions to the binding energy of ID CIFMs. Finally, 
our main results are summarized in Sec. [VI] 



II. TWO-CHANNEL HAMILTONIAN 

For a two-component Fermi gas with atomic mass m, 
the two-channel effective Hamiltonian that we consider 
includes the following single-particle Hamiltonian and in- 
teraction Hamiltonian [2ll |22j , 



Ho = I ( r 
d 3 r¥ (r) 
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and 



H in t = U J d 3 r^\ (r) i,\ (r) ^ (r) ^ (r) 

+ g f d 3 v {*t ( r ) ^ ( r ) ^ ( r ) + J . ( 2 ) 



Here, t/j a (r) (a I) and ^ (r) are the field operators 
of fermionic atoms in the open (atomic) channel and 
of bosonic molecules in the closed (molecular) channel, 
respectively. V ex t (r) is the external trapping potential, 
which could be described by a harmonic trap as a good 
approximation. In order to realize a 2D system, a tight 
axial confinement is applied experimentally, while the ra- 
dial (transverse) confinement is much weaker. Thus, the 
radial motion of atoms is approximatively free. Theo- 
retically, we can only consider a tight axial confinement 
V e xt ( z ) = muj 2 z 2 /2 along the z-direction. In a like man- 
ner, for a ID system, the external trapping potential can 
be chosen as V ex t (p) = mw 2 _p 2 /2, which is in x — y plane 



with p 



-y . A (B) is the detuning between atomic 



and molecular channels, which could be tuned by a mag- 
netic field near a Feshbach resonance. U describes the 
bare interaction between atoms with different spins in 
the atomic channel, and g denotes the coupling strength 
between the atomic and molecular channels. We have 
used the constants U and g to describe the coupling pro- 
cesses both for the bare interatomic interaction and the 
formation of molecules, which is analogous to those used 
in the pseudopotential theory. Consequently, a diver- 
gence in high energy will appear. Thus, a regularization 
must be introduced as we shall mention later. 



III. ANSATZ OF THE TWO-BODY 
WAVEFUNCTION 



In order to calculate the binding energy of CIFMs, let 
us solve a two-body problem. In this paper, we are only 
interested in the tight- confinement limit, htu\\ t ± S> fc^T , 
where ks is the Boltzmann constant and T is the temper- 
ature. For an ultracold Fermi gas with a pretty low tem- 
perature, almost all the atoms occupy the lowest-energy 
(ground) state of the confined direction and can not be 
excited to the higher energy levels during collisions. Due 
to the harmonic confinement, the COM motion is com- 
pletely decoupled from the relative motion of two atoms. 
Therefore, the COM of the two atoms will always stay in 
the ground state of the confined direction. In addition, 
there is always a global translation invariance in the un- 
trapped directions, which means the COM momentum 
K in these directions is a good quantum number. Thus, 
an ansatz of the two-body wavefunction could take the 
following form, 



j d 3 re lKp <p {z) & (r) 
J d 3 v l d 3 T 2 ^^+P^I 2 ^ 



z\ + z 2 



xQ 2 (r 1 -ra)^{(r 1 )^J(r a ) |0> (3) 
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for a 2D system, and, 

|^) = [/ d 3 re iKz <p 00 (p)&(v) + 
d 3 rid 3 r 2 e iK(zi+22)/2 0oo 



function © is valid only for e < 0. Obviously, when 
r — » 0, the Green's function diverges as, 



P1+P2 



xQiCri-^)^^)^^)] |0) (4) 

for a ID system. Here |0) stands for the vacuum state, 
and Q^) (r) is the relative wavefunction of the two atoms 
in a "D"-dimensional system, which is spatially antisym- 
metric, e.g., Q 2 (-r) = -Q 2 (r). O o (p) and 4>o(z) are 
the ground states of ID and 2D harmonic oscillators, re- 
spectively. 



IV. SOLUTION OF TWO-BODY PROBLEM 

A. In a 2D system 

Let us first consider an axially confined Fermi gas, 
where atoms can move freely in the x — y plane. By 
acting the two-channel Hamiltonian (pj and ([2]) onto the 
two-body wavefunction ansatz the Schordinger equa- 
tion H \(p 2 ) = £ \f2) is equivalent to the following set of 
equations (see the Appendix for details), 
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which results in an energy divergence. This is what we 
have anticipated since we use a pseudopotential method 
to describe the coupling process both for the interatomic 
interactions and for the formation of molecules. 

In order to eliminate this singularity, the Green's func- 
tion must be regularized as follows, 



lim £ £ (2) (r, 0) = lim 

1 — >0 r— >0 
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Then, by substituting Eqs.flTT]) and ([12]) into Eq.©, 
Q 2 (r)+[J7Q 2 (r) + g] S (r) = EQ 2 (r) ,Q 2 (0) can be easily solved, 



E = A (B) + gQ 2 (0) 



(3) 



(6) 



Q 2 (0) 



2ir 3/2 d\huj\\ T 2 X (e) + £/ 
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where the COM part has been separated out. r = ri — r 2 
is the relative coordinate, [i = m/2 is the reduced mass, 
and E = e — H 2 K 2 /4m — fiwy/2 is the relative energy. 
From Eq.([5]), the bound-state solution can be written as, 



Q 2 (r)=-[UQ 2 (0)+g] G% ] (r,0), 



(7) 



(2) 1 

where G E (r, r') is the Green's function that satisfies, 



2i f7 2 + ^p 2 -E)Gf(r.r')= t )ir r'i. (S) 



After some straightforward algebra [L 
(2) 

function G E (r, 0) is easily obtained, 



the Green's 
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Consequently, by combining Eqs.© and (fTB"]) . the bind- 
ing energy Eb of 2D CIFMs should satisfy the following 
self-consistent equation, 



E B = A (S)- 5 2 27r 3/2 d|fto;|| J"^ 1 (E B /Hu\\ - 1/2) + 17 

(14) 

By solving this self-consistent equation, the binding en- 
ergy can be obtained. 



B. In a ID system 

For a ID system, atoms can only move freely along the 
axial direction, while the radial motion is frozen. Oper- 
ating on \ipi) , Eq.((4]), with the Hamiltonian JTJ) and |[2j), 
the Schordinger equation is deduced to the following set 
of equations, 



-|-V 2 + l^lp 2 \ Q x (r)+[UQx (r) + g] S (r) = EQx (r) 



tvT 



(9) 



2/t 



where dii = * t[ /h/ fMjJn is the harmonic length, and e = 
Ejhjj\\ — 1/2 . This integral representation of the Green's 



E — A (B) + gQx (0) 



(15) 
(16) 
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where the relative energy is E = e — h 2 K 2 /4m — haj±. 
After the similar derivation as that for the 2D self- 
consistent equation of the binding energy Eb can easily 
obtained, 

E B = A (B)-g 2 2vr 3/2 difiw ± J"f 1 (E B /fru± - 1) + U 

(17) 

where dj_ = \fh[Jujj±_ and 



exp(ei/2) 1 
Vi(l-e-*) ~ t^ 2 



(18) 



Then the binding energy of the ID CIFMs can be solved 
using the self-consistent equation (fl7|). 
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V. RESULTS AND DISCUSSION 

A. At the resonance 

The single-channel theory is used to describe the uni- 
versality for a strongl y- in teracting Fermi gas near a broad 
Feshbach resonance [23l l24l |. in which all other length 
scales become irrelevant except the average interatomic 
distance. However, there is always an intrinsic length 
r* for all Feshbach resonances, which is related to the 
resonance width. It can be defined as (l9| . 



(19) 



2mab g AiiAB 



where ab g is the background scattering length in the open 
channel, AB is the width of the resonance, and A/i is 
the magnetic moment difference of the molecular state 
with respect to the threshold of two free atoms. For a 
broad resonance, this intrinsic length is quite small due 
to the large resonance width, thus it is reasonable that 
the single-channel theory gives a good description of Fes- 
hbach resonances. However, when the resonance width 
becomes narrow, the intrinsic length r* can not be ig- 
nored anymore, and will obviously affect the properties 
of the resonance. Consequently, the prediction of the 
two-channel theory, which is physically more realistic, is 
expected to deviate from that of the single-channel the- 
ory. In order to demonstrate this deviation, let us focus 
on the situation right at the resonance, where the bare 
interatomic interaction U could be reasonably neglected, 
as we shall discuss in greater detail later. The channel 
coupling strength g is given by (l8l [25| , 



g = h 



4na bg ABAfi _ H 2 



(20) 



Combining with Eq. tfMf . the binding energy of 2D 
CIFMs predicted by the two-channel theory at the reso- 
nance can be written as, 



'II 
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Figure 1: (Color online) The binding energy of 2D CIFMs 
as a function of d\\/r* right at the Feshbach resonance. The 
red solid circles with error bar are taken from the 2D 40 K 
experiment near the Feshbach resonance at Bo = 202. 1G. 
The two-channel predictions for the binding energy of two 
relatively narrow resonances for 2S Na at Bo = 907G (blue 
circle) and 87 Kb at Bo = 1007. 4G (wine triangle) are also 
presented. Here, the trap frequency cjy is chosen to be 2n x 75 
kHz as in Ref. 



where eb = Eb — . 

The binding energy e b / ip^\\) of 2D CIFMs as a func- 
tion of the intrinsic length dii/r* is shown in FigfJJ At 
broad resonance limit, i.e., dii/r* — > oo , the two-channel 
prediction approaches to the result of the sin gle- channel 
theory, which is predicted to be — 0.244fiw|| [26|. How- 
ever, as the resonance becomes narrow, the two-channel 
prediction will obviously deviate from the single-channel 
result, and the binding energy of 2D CIFMs is strongly 
dependent on the intrinsic length or the resonance width. 
In the recent 2D 40 K experiment [17| . the binding energy 
of 2D CIFMs is measured across a broad Feshbach reso- 
nance at Bq = 202. 1G with dn/r* ~ 110. The observed 
binding energy near the resonance is indicated in Figfflby 
the red solid circles with error bar. Since there is no ex- 
perimental data exactly on the resonance, two points near 
the resonance are chosen, and the binding energy right 
at the resonance should lie between these two points. 

We also predict the binding energy of 2D CIFMs for 
two relatively narrow resonances of 23 Na at Bq = 907G 
and 87 Rb at B = 1007. 4G, which are presented by the 
blue circle and wine triangle in FigO] respectively. 

In a like manner, the binding energy of ID CIFMs at 
the resonance is given by, 



( e B 



- 1 



(22) 



where eb = E — fhu)± . The binding energy of ID CIFMs 
as a function of the intrinsic length is plotted in FigJ5J 
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Figure 2: (Color online) The same as Fig[TJbut for ID CIFMs. 
The red solid circles with error bar are taken from the ID 40 K 
experiment showing the two data points nearest to the 
Feshbach resonance. Here, the trap frequency lj± is chosen 
to be 2tt x 69 kHz as in Ref. @|. 



B. 2D CIFMs across the Feshbach resonance 
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Figure 3: (Color online) The binding energy of 2D CIFMs 
across the Feshbach resonance at Bo = 202. 1G, comparing 
with the experimental data (solid circles) fl7l |. The param- 
eters for this resonance are AB — 7.8G, Afj, — 1.68/is, 
aig = 174as, where fJB and as are the Bohr magneton and 
radius, respectively. The trap frequency in the strongly con- 
fined direction is u\\ = 2n x 75 kHz. 



In the following, let us investigate the binding energy 
of 2D CIFMs across Feshbach resonances. By using Eqs. 
(fT2"]) and HU), the binding energy predicted by the two- 
channel theory satisfies, 



e B + -^H = A{B)~g 2 



1 njjj\ 



U 



(23) 

We recall that the binding energy given by the single- 
channel theory is determined by the following equation 



d 
a3D 



Tiw\\ 



(24) 



where a 3 D is the effective 3D s-wave scattering length. 

In a magnetic Feshbach resonance, the detuning 
between the open and closed channels, A (B) = 
Afj, (B — Bo) , varies as a function of the magnetic field B 
and depends on the moment difference Afj, between the 
open and closed channels. The bare interatomic inter- 
action U and the channel-coupling constant g a re g iven 
by U = 4nh 2 a bg /m and Eq-PJ]), respectively ^ HI- 
The effective 3D scattering length a 3 u is adjusted by the 
magnetic field B according to, 



0-3D (B) = a bg 1 



AB 

B-B 



(25) 



The two-channel prediction of the binding energy of 2D 
CIFMs across a broad Feshbach resonance as well as that 
of the single-channel theory is illustrated in Fig EH com- 
paring with the recent 2D 40 K experiment [17|. It can 
be seen clearly that the two-channel theory agrees with 



the single-channel result for such a broad resonance, and 
both agree well with the experiment at the resonance. 
However, away from the resonance, the observed bind- 
ing energy deviates from both two-channel and single- 
channel predictions. 



C. ID CIFMs across the Feshbach resonance 

We now turn to the binding energy of ID CIFMs. Ac- 
cording to Eqs. pT|) and (fT5|) . the two-channel prediction 
for the binding energy of ID CIFMs is given by, 



€ B + huj ± = A (B) — g 2 



2TT 3/2 d 3 L h0J_ L X 



u 



(26) 



The two-channel result is plotted as well as the single- 
channel prediction in Fig. [4] We find that the results of 
the two-channel and single-channel theories are almost 
the same, and both agrees with the experiment [6j at 
the resonance. This agreement is anticipated since the 
resonance used in Q is also very broad, and the two- 
channel theory should approach the single-channel result, 
as we already seen in Fig. [21 

In Ref. , Dickerscheid and Stoof suggested to ignore 
the bare interatomic interaction U , by assuming that this 
background interaction is quite small comparing to the 
channel coupling near the resonance. Under this approx- 
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Figure 4: (Color online) The binding energy of ID CIFMs 
across the Feshbach resonance at Bo = 202. 1G, comparing 
with the experiment data (solid circles) Q. The blue dash- 
dotted curve is the result of Ref . P3 , obtained by ignoring the 
bare interatomic interaction in the open channel, i.e., U — 0. 
Here, the trap frequency in the strongly confined direction is 
oj ± = 2tt x 69 kHz Q. 



imation, Eq. ([26)1 yields 
e B + hu± ~ A (B) — 



27r 3 / 2 d3 hw 




(27) 



Here, we de- 
H//mj± , which 



which recovers the results in Ref. 
fine a different harmonic length dj_ - 
is V2 times larger than that in Ref. [18|, and J-\ (e) 
V^C (1/2, — e) , where C{s,q) is the Hurwitz Zeta func- 
tion. The results of Eq. ([27|) is plotted in Fig. H by a 
blue dash-dotted curve. It seems that the two-channel 
theory provides an improved agreement with the exper- 
imental data, compared with the single-channel theory. 
Here, however, we find that if the bare interatomic in- 
teraction is correctly included, the result of the full two- 
channel theory returns back to that of the single-channel 
calculation, as expected for a broad Feshbach resonance. 

Now, it is interesting to discuss the condition under 
which the bare interatomic interaction U could be ne- 
glected. From Eq. (|26[) . we find whether the bare inter- 
atomic interaction could be neglected lies on the competi- 
tion between 2tt 3 '' ' 2 d^tkJxJ 7 ^ [es/ (hu)±)] and U, rather 
than the ratio of g/U. Thus, the condition can be taken 



» U, 



or, 



( e B 



a bg 



(28) 



(29) 



The function T\ [cb/ (Huj±)] as a function of the mag- 
netic field is presented in Fig. [S] Obviously, the condi- 
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Figure 5: (Color online) The function T\ [es / (hoj±)] , in units 
of y/Trd±/ab g , varies as a function of the magnetic field for 
a fD system across the resonance Bo = 202. fG. The dashed 
line indicates the location of the Feshbach resonance. 



tion (|29j) can be fulfilled well very close to the resonance, 
which means ignoring the bare interatomic interaction 
is indeed a good approximation at the resonance limit. 
However, away from the resonance, we see that the bare 
interatomic interaction will bring a sizable correction to 
the binding energy of the molecules, which makes the 
approximation results deviate from those of the full two- 
channel calculations, as we have already seen in Fig. 2) 



VI. CONCLUSIONS 

In summary, we have presented a two-channel theory 
for calculating the binding energy of confinement-induced 
Feshbach molecules in 2D and ID ultracold atomic sys- 
tems. We have found that the two-channel results will de- 
viate from the single-channel predictions as the width of 
Feshbach resonances decreases, although these two the- 
ories give the same results in the limit of broad reso- 
nances. Compared with the recent experiments, where 
a broad magnetic Feshbach resonance is used to control 
the interatomic interactions, both the two-channel and 
single-channel theories agree well with the data near res- 
onance. 

For ID confinement-induced Feshbach molecules, a full 
two-channel calculation has been performed. Compared 
with the previous work [18j, in which the background 
interatomic interaction in the open channel is neglected, 
we have found that the bare interatomic interaction will 
yield a sizable correction to the binding energy away from 
the resonance. Thus, in contrast to the difference shown 
in Ref. [l8[ , the two-channel result agrees very well with 
the single-channel prediction, as we may anticipate for a 
broad Feshbach resonance. 
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Appendix : The Schordinger equation for a 2d system 



In this appendix, we present the detailed derivation of how the 2D Schordinger equation % \ip 2 ) — e \tp 2 ) is reduced 
to the coupled equations ([5]) and ([6]). By using the commutation and anti-commutation relations of the Bose and 
Fermi field operators. 



and 



$ (r) , (r') =6(r-r'), 
^(r),V^(r')l = <5(r-r')<W, 



(30) 
(31) 



where =p represent the commutation and anti-commutation relations, respectively, we obtain, 
Ufa) = yd 3 rid 3 r 2 ^(n)^j(r 2 ) 



h 2 h 2 

^V 2 ri + V ext (n) - — V 2 2 + V ext (r 2 ) + UiS (n - r 2 ) 



3l K(p 1+P2 )/ 20o (^±^l\ Q 2 (ri - r2 ) + gS ( ri - r 2 ) e *(*+ft)/% (^4^) } 



|0) 



d 3 r& (r) 



*h 2 K 2 1. 



4m 



-hu + A{B)+gQ 2 (0) 



e iK '^o (z) |0) 



(32) 



Here, we have introduced pseudopotentials to describe the bare interatomic interaction and the channel coupling. 
Then, by comparing the corresponding terms in % \ip 2 ) and e \f 2 ), we arrive at, 



h 2 h 2 

■TT-Vti + V *** ( ri ) " ^T- V r 2 + V** (r 2 ) + U5 (n - r 2 ) 



3 ,K(p 1+P2 )/ 20o ^fl + i2^ g2(ri _ r2) 



+ <?<5(r 1 -r 2 )e iK ^+^)/ 2 o 



z\ + z 2 



ee4 K(p 1+ p 2 )/2 0o + } Q2 (ri _ r2) ; (33) 



and 



*h 2 K 2 1 



fiwii + A (B) + gQ 2 (0) 



iKp 



0o («) = ee lK -"0 o (z) 



4m 2 

If we seperate the COM motion from the relative motion, Eqs. ([3"3"|) and (|3~4"|) yield, 

h 2 1 

' t -72 , 1 ,,, ,2 2 



and 



— V 2 + -^\z 2 ) Q 2 (r) + [UQ 2 (r) + <?] <5 (r) = EQ 2 (r) 



(34) 



(35) 



E = A (B) + gQ 2 (0) 



(36) 

respectively, where r = — r 2 , /i = m/2 is the reduced mass, and E = e — h 2 K 2 /4m — Hu)\\/2 is the relative energy. 
In a like manner, we can easily obtain Eqs. flTSj) and (flo| for ID systems. 
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